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AERONAUTIC  SYMBOLS 
1,  FUNDAMENTAL  AND  DERIVED  UNITS 


Symbol 

.  Metric 

English 

Unit 

Abbrevia¬ 

tion 

Unit 

Abbreviation 

Length _ 

Time--  --  -_ 
Force _ 

1 

t 

F 

meter-  -_  _  _ 

second _  _  - _ 

weight  of  1  kilogram.-  .. 

m 

s 

kg 

foot  (or  mile)  .  _  -  -  - 

second  (or  hour)  _ 

weight  of  1  pound _ 

ft  (or  mi) 
sec  (or  hr) 
lb 

Power _ 

P 

V 

horsepower  (metric) _ 

horsepower 

hp 

mph 

fps 

Speed - 

/kilometers  per  hour _ 

\meters  per  second  _  -  - 

kph 

mps 

miles  per  hour _ 

feet  per  second-- 

W 

g 

m 

I 

s. 

G 

b 

c 

A 

V 

g 

L 

D 

D, 

Dt 

D, 

C 


2.  GENERAL  SYMBOLS 


Weight = mgr 

Standard  acceleration  of  gravity= 9.80665  in/s^ 
or  32.1740  ft/sec^ 

W 

Mass=— 

g 

Moment  of  inertia=mP.  (Indicate  axis  of 
radius  of  gyration  k  by  proper  subscript.) 
Coefficient  of  viscosity 


V  Kinematic  viscosity 

p  Density  (mass  per  unit  volume) 

Standard  density  of  dry  air,  0.12497  kg-m~^-s^  at  15°  C 
and  760  mm;  or  0.002378  Ib-ff*  sec^ 

Specific  weight  of  “standard”  air,  1.2255  kg/m®  or 
0.07651  Ib/cu  ft 


3.  AERODYNAMIC  SYMBOLS 


Area 

Area  of  wing 
Gap 
Span 
Chord 

Aspect  ratio,  ^ 

True  air  speed 
Dynamic  pressure,  ^  pK® 


Lift,  absolute  coefficient 


D 


Drag,  absolute  coefficient 
Profile  drag,  absolute  coefficient 


D, 


Induced  drag,  absolute  coefficient 


Parasite  drag,  absolute  coefficient  Goj,—  ^ 


Cross-wind  force,  absolute  coefficient  Cc 


C 

"qS 


in  Angle  of  setting  of  wings  (relative  to  thrust  line) 
it  Angle  of  stabilizer  setting  (relative  to  thrust 
line) 

Q  Resultant  moment 

Q  Resultant  angular  velocity 

R  Re3molds  number,  p  ^  where  f  is  a  linear  dimen¬ 
sion  (e.g.,  for  an  airfoil  of  1.0  ft  chord,  100 
mph,  standard  pressure  at  15°  C,  the  corre¬ 
sponding  Reynolds  number  is  935,400;  or  for 
an  airfoil  of  1.0  m  chord,  100  mps,  the  corre¬ 
sponding  Reynolds  number  is  6,865,000) 
a  Angle  of  attack 

i  Angle  of  down  wash 

ao  Angle  of  attack,  infinite  aspect  ratio 

a,;  Angle  of  attack,  induced 

a.  Angle  of  attack,  absolute  (measured  from  zero- 

lift  position) 

y  Flight-path  angle 
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ANALYTICAL  DETERMINATION  OF  COUPLED  BENDING-TORSION  VIBRATIONS 
OF  CANTILEVER  BEAMS  BY  MEANS  OF  STATION  FUNCTIONS' 

By  Alexandeb  Mendelson  and  Sblwyn  Gendleb 


SUMMARY 

A  method  based  on  the  concept  of  Station  Functions  is  pre¬ 
sented  for  calculatmg  the  modes  and  the  frequencies  of  non- 
uniform  cantilever  beams  vibrating  in  torsion,  bending,  and 
coupled  bending-torsion  motion.  The  method  combines  some 
of  the  advantages  of  the  Rayleigh-Ritz  and  Stodola  methods,  in 
that  a  continuous  loading  function  for  the  beam  is  used,  with 
the  advantages  of  the  influence-coefficient  method,  in  that  the 
continuous  loading  function  is  obtained  in  terms  of  the  dis¬ 
placements  at  a  finite  number  of  stations  along  the  beam. 

The  Station  Functions  were  derived  for  a  number  of  stations 
ranging,  from  one  to  eight.  The  deflections  were  obtained  in 
terms  of  the  physical  properties  of  the  beam  and  Station  Num¬ 
bers,  which  are  general  in  nature  and  which  have  been  tabulated 
for  easy  reference.  Examples  were  worked  out  in  detail;  com¬ 
parisons  were  made  with  exact  theoretical  results.  For  a  uni¬ 
form  cantilever  beam  with  n  stations,  the  first  n  modes  and 
frequencies  were  in  good  agreement  with  the  theoretically  exact 
values.  The  effiect  of  coupling  between  bending  and  torsion 
was  shown  to  reduce  the  first  natural  frequency  to  a  value  below 
that  which  it  would  have  if  there  were  no  coupling. 

INTRODUCTION 

The  failure  of  turbine  and  compressor  blades  due  to  vibra¬ 
tions  has  led  to  an  increased  interest  in  the  study  of  the 
vibrations  of  these  blades  and  in  the  determination  of  the 
natural  modes  and  frequencies.  In  such  theoretical  studies, 
it  is  usually  assumed  that  the  compressor  or  turbine  blade 
acts  as  a  cantilever  beam.  The  calculation  of  the  imcoupled 
modes  of  arbitrarily  shaped  cantilever  beams  has  been  ex¬ 
tensively  investigated  (references  1  to  4),  but  little  work  has 
as  yet  been  done  on  calculating  the  coupled  modes  of  such 
beams.  If  the  geometry  of  the  beam  is  such  that  coupling 
exists,  the  coupled  modes  are  the  actual  vibrational  modes 
that  must  be  calculated. 

Four  general  methods  are  currently  in  use  for  calculating 
uncoupled  modes  and  frequencies  of  nonuniform  beams. 
These  methods  are  the  Rayleigh-Ritz  or  energy  method 
(reference  1),  the  Stodola  method  (references  5  and  6),  the 
influence-coefficient  method  (references  4  and  7),  and  the 
integral-equation  method  (references  8  and  9).  For  each  of 
these  methods,  computational  work  can  usually  be  carried  out 
in  several  ways.  For  example,  by  the  use  of  influence  co¬ 
efficients  the  modes  and  frequencies  can  be  detenftined  by 

>  Supersedes  NACA  TN  2185,  “Analytical  Determination  of  Coupled  Bending-Torsion  V 
Selwyn  Gendler,  1950. 


Mykelstad’s  iteration  procedure  (reference  7)  or  by  matrix 
methods  (reference  4). 

Any  one  of  these  methods  can  be  extended  to  the  calcula¬ 
tion  of  coupled  bending-torsion  modes.  The  Rayleigh-Ritz 
method  usually  requires  that  the  uncoupled  modes  be  deter¬ 
mined  before  the  coupled  modes  can  be  computed.  In  apply¬ 
ing  either  the  Rayleigh-Ritz  or  the  Stodola  method,  great 
difficulty  is  encoimtered  in  accurately  determining  the  higher 
modes,  because  the  lower  modes  must  first  be  “swept  out” 
by  the  use  of  exact  orthogonality  conditions  (reference  10); 
the  process  will  otherwise  always  converge  back  to  the 
lowest  mode.  The  same  difficulties  are  encountered  in  the 
integral-equation  method. 

The  influence-coefficient  method  reduces  the  problem  to 
one  having  a  finite  number  of  degrees  of  freedom.  The  beam 
is  divided  into  n  intervals  and  a  concentrated  loading  is  as¬ 
sumed  at  the  center  of  gravity  of  each  interval.  The  solution 
of  the  resultant  determinantal  equation  gives  the  first  n 
modes.  The  accuracy  of  the  higher  modes  is,  however,  very 
poor;  only  the  first  third  of  the  modes  and  the  first  half  of  the 
frequencies  are  obtained  within  the  usual  engineering  accu¬ 
racy.  Carrying  along  so  many  useless  modes  greatly  in¬ 
creases  the  labor  involved. 

A  straightforward  accurate  method  for  determining  the 
coupled  bending-torsion  modes  and  the  frequencies  of  non- 
uniform  cantilever  beams,  together  with  applications  of  this 
method,  was  developed  at  the  NACA  Lewis  laboratory  dur¬ 
ing  1949  and  is  presented  herein.  This  method  is  based  on 
the  use  of  Station  Functions  as  first  discussed  in  refer¬ 
ence  11.  Incorporated  in  the  method  are  the  advantages  of 
the  continuous-function  deflections  of  the  Rayleigh-Ritz 
and  Stodola  methods  together  with  the  advantages  of  the 
finite  number  of  degrees  of  freedom  of  the  influence-coefficient 
method.  When  the  method  is  applied  to  a  uniform  beam, 
the  first  n  roots  of  the  resultant  determinantal  equation  are 
amply  accurate  for  engineering  purposes. 

The  final  determinantal  equation  is  solved  herein  by 
matrix-iteration  methods  (reference  4).  Any  other  con¬ 
venient  method  may,  however,  be  used  and  no  knowledge 
of  matrix  algebra  is  needed  to  carry  out  the  calculations  by 
the  matrix  method.  The  work  can  be  done  by  an  inexperi¬ 
enced  computer,  as  the  only  operations  necessary  for  determ¬ 
ining  each  mode  are  cumulative  mifitiplication  and  division. 
In  addition,  for  the  case  in  wliich  the  coupling  coefficient 
remains  constant  along  the  beam,  a  simple  quadratic 
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formula  and  a  series  of  curves  are  presented  for  determining 
the  first  coupled  mode  in  terms  of  the  uncoupled  modes. 
Examples  are  developed  in  detail  and  comparisons  with 
exact  theoretical  results  are  included. 

THEORY 

In  the  usual  influence-coefficient  methods  for  solving 
dynamical  problems,  a  continuous  body  having  an  infinite 
munber  of  degrees  of  freedom  is  replaced  by  a  body  having  a 
finite  number  of  degrees  of  freedom.  Two  principal  assump¬ 
tions  are  then  made  that  introduce  inaccuracies  into  the 
solutions,  particularly  in  the  higher  modes:  (1)  The  resultant 
of  the  inertia  loads  of  all  the  infinitesimal  masses  in  a  finite 
interval  passes  through  the  center  of  gravity  of  that  interval; 
and  (2)  a  concentrated  load  that  is  the  resultant  of  a  dis¬ 
tributed  load  produces  the  same  deflection  as  the  distributed 
load.  An  attempt  has  been  made  to  reduce  the  error  due  to 
the  second  of  these  assumptions  by  the  use  of  weighting 
matrices  (reference  12).  Although  the  accuracy  is  thereby 
increased,  the  effect  of  the  first  assumption  is  still  great 
enough  to  introduce  serious  errors  (reference  11). 

In  order  to  eliminate  these  assumptions,  Rauscher  (ref¬ 
erence  11)  introduced  the  concept  of  Station  Functions. 
Instead  of  assuming  the  inertia  loads  to  be  concentrated  at 
the  centers  of  gravity  of  the  intervals,  the  inertia  loads 
and,  consequently,  the  deflections  are  assumed  to  be  con¬ 
tinuous  functions  along  the  beam.  The  values  of  these 
continuous  deflection  functions  at  the  reference  stations  must 
equal  the  deflections  of  the  reference  stations.  The  loading 
on  the  beam  is  therefore  a  continuous  fimction  of  the  de¬ 
flections  of  the  reference  stations.  Inasmuch  as  the  deflec¬ 
tions  of  the  reference  stations  can  be  computed  from  the 
loading  on  the  beam,  which  in  turn  is  available  from  the 
deflections,  the  deflections  are  therefore  obtained  as  functions 
of  themselves.  This  procedure  gives  n  homogeneous  equa¬ 
tions  in  the  n  deflections  of  the  reference  stations.  The 
resultant  determinantal  equation  has  n  roots  for  the  fre¬ 
quency;  it  will  be  shown  that  for  a  uniform  beam  all  these 
roots  are  sufficiently  accurate  for  engineering  purposes  if 
the  deflection  functions  are  properly  chosen.  (For  coupled 
bending-torsion  vibrations,  2n  homogeneous  equations  and 
2n  roots  are  obtained  for  n  stations.) 

The  deflection  functions  used  must  satisfy  the  boundary 
conditions  of  the  problem  and  also  the  condition  that,  at 
any  reference  station,  the  value  of  the  function  must  equal 
the  deflection  of  the  reference  station.  Although  it  is  always 
possible  to  find  directly  a  single  fimction  that  will  satisfy 
these  conditions,  it  is  more  convenient  to  obtain  different 
component  functions  at  each  station  and  to  add  all  these 
component  functions  together  to  give  the  complete  deflec¬ 
tion  function.  Rauscher  (reference  11)  calls  these  compon¬ 
ent  deflection  functions  Station  Functions.  For  example, 
the  complete  torsional  deflection  function  for  the  beam  will 
have  the  following  form : 


jj{z)  Station  Function  in  torsion  associated  with  station 
(All  symbols  are  defined  in  appendix  A.) 

Each  Station  Function  must  satisfy  the  boundary  condi¬ 
tions  of  the  problem  and  the  following  additional  conditions : 
(1)  At  the  reference  station  with  which  it  is  associated,  the 
Station  Function  equals  the  deflection  of  that  reference  sta¬ 
tion;  and  (2)  at  all  other  reference  stations,  the  Station 
Function  equals  zero.  The  sum  of  all  these  Station  Func¬ 
tions  will  then  give  the  complete  deflection  function  for  the 
beam.  The  Station  Functions  and  corresponding  loading 
functions  are  derived  in  appendix  B  for  torsional  vibrations. 
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Figure  1.— Cantilever  beam  with  n  stations. 


bending  vibrations,  and  coupled  bending-torsion  vibrations  of 
an  arbitrary  cantilever  beam. 

Torsional  vibrations. — It  is  shown  in  appendix  B  that  the 
torsional  deflections  of  the  reference  stations  for  a  beam 
divided  into  n  intervals  of  length  5,  as  shown  in  figure  1, 
are  given  by  the  following  system  of  equations : 

Tr.  ™ 

(1) 

Uo  j=i 

where 

aij^±,)rlhNj,-{k-l)hM,,+  ±,  /.A/J  (2) 

4  =  1  Uj  L  r=4+l  J 

i  and  1, 2,  .  .  .  n 


CO  frequency  of  vibration 

S  length  of  interval 

Iq  mass  moment  of  inertia  per  unit  length  about  elastic 
axis  at  root  section 

If,  ratio  of  average  mass  moment  of  inertia  per  unit  length 
of  interval  to  mass  moment  of  inertia  per  unit  length 
at  root  section 

Co  torsional  stiffness  of  root  section 

Cf  ratio  of  average  torsional  stiffness  of  A:®  interval  to  tor¬ 
sional  stiffness  at  root  section 
The  Station  Numbers  Nff  and  Mjf  are  functions  only  of  the 

integers  k,  j,  and  n  and  are  defined  as 


(3) 


where  Jj{z)  represents  the  Station  Functions  derived  in 
appendix  B  and  is  given  by 


(03) 

J=1 

where  *•  .  ’ 

z  dimensionless  distance  along  beam 

e{z)  torsional  deflection  at  distance  z  from  root 
Bj  torsional  deflection  at  station 


ji{z)^aiiZ-\-a2jZ‘^-\- .  .  .  (4) 

The  coefficients  a,,  are  determined  in  appendix  B  by 
satisfying  the  conditions  on  the  Station  Functions.  The 
integrals  in  equations  (3)  are  thus  seen  to  be  integrals  of 
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simple  polynomials  and  the  limits  of  integration  are  integers. 
The  Station  Numbers  and  Mji,  are  therefore  rational 
numbers,  functions  only  of  the  integers  n,  k,  and  j.  These 
numbers  have  been  evaluated  and  are  listed  in  tables  I  to 
VIII. 

If  the  physical  properties  of  the  beam  under  consideration 
are  known  for  each  of  the  n  intervals,  and  /*  will  be 
known.  The  Station  Numbers  Nj^  and  Mjt  can  be  obtained 
from  tables  I  to  VIII.  From  equation  (2),  aij  can  then  be 
easily  calculated. 

Equation  (1)  actually  represents  n  homogeneous  equations 

1  O 

in  the  n  unknown  deflections  0,-.  With— 2^^=X,  these  equa- 

Cl) 


(5) 


tions  can  be  written  as  follows: 

(an  —  X)  ^i  +  ai2^2  +  «i3^3  + 
a21  ^1  +  (a22  —  X)  ^2  +  «23^3  + 

«31^1  +  «32®2+(“33  —  X)03  + 

.  .  .  ~\~  OCinOn  —  0 

■  •  •  -\-OC2n0n=O 

■  ■  ■  +a3n^«  =  0 

«Ml®l  +  n:n2^2  +  aM3^3+  •  •  • 

+  (anK  —  X)0„  =  O 

For  a  nontrivial  solution,  the  determinant  of  the  coefficients 
must  vanish  and  the  characteristic  equation  becomes 


an  — 

X  aj2 

“13 

•  •  “in 

a2i 

a22— 

CO 

.•  •  a2n 

asi 

“32 

“33  ^  • 

•  .  “3n 

=  0 

(6) 

a«i 

an2 

“«3  •  • 

•  ^nn  ^ 

[X/— [a<y]l==0  (6a) 

where  I  is  the  identity  matrix,  and  [ay]  is  the  dynamical 
matrix. 

Equation  (6)  can  be  solved  for  the  n  values  of  X  by  any 
method  available.  The  method  used  herein  was  to  obtain 
the  values  of  X  as  the  latent  roots  of  the  matrix  [ay],  which 
is  actually  the  dynamical  matrix  for  the  problem.  The  mode 
shapes  are  obtained  at  the  same  time. 

Bending  vibrations. — The  bending  deflections  for  the  beam 
shown  in  figure  1  are  given  by  the  following  system  of  equa¬ 
tions  (appendix  B) : 


2/j  =  a)^5^^  S  ^uV) 

■Oo  j=l 


(7) 


where 


+ 


[ 


F-(it-l)*  i2k-l)i 


0  ^'"0 


(8) 


i  and  j—l,2,  .  .  .  n 

mo  mass  per  imit  length  of  beam  at  root  section 
Mk  ratio  of  average  mass  per  imit  length  of  k*^  interval  to 
mass  per  unit  length  at  root  section  , . 

Bq  bending  stiffness  at  root  section 

Bk  ratio  of  average  bending  stiffness  of  k*^  interval  to  bend¬ 
ing  stiffness  at  root  section 


The  Station  Numbers  M'j^,  N'jjc,  P'jh,  and  Q' jk  are  func¬ 
tions  only  of  the  integers  k,  j,  and  n  and  are  defined  by 

=  1)2+1  gi{z)d2 


+l  (^  —  1)^1  9j{z)di 


(9) 


M'ik=  f  gAz)dz 

Jk-l 

N'jk=  f  zgj{z)dz 
Jk-l 

The  Station  Functions  gj{z)  are  derived  in  appendix  B  and 
are  given  by 

gj{z)==b2jZ^+b3jZ^+bijZ^+  .  .  .  +6(„+3)j2'“+®>  (10) 

The  integrals  in  equations  (9)  are  thus  seen  to  be  integrals 
of  simple  polynomials.  The  Station  Numbers  M'jic,  N'jk, 
P' jk,  and  Q! jk  are  rational  numbers,  functions  only  of  the 
integers  j,  k,  and  n.  These  numbers  have  been  evaluated 
and  are  listed  in  tables  I  to  VIII. 

If  the  physical  properties  of  the  beam  are  known  for  each 
of  the  n  intervals,  and  Bk  will  he  known.  The  Station 
Numbers  M'jk,  N'jk,  P'jk,  and  Q'jk  are  obtained  from  tables  I 
to  VIII;  |3y  can  then  easily  be  calculated  by  equation  (8). 
The  determinantal  equation  is: 


fln-X 

fll2 

fll3  . 

■  flin 

fl21 

fl22  — 

X  fl23 

■  fl2» 

fl31 

fl32 

flss  X  . 

•  flsn 

=0 

(11) 

flnl 

fln2 

flnS 

•  ^nn  ^ 

lx/ — [fly]!  —  0 


(11a) 


where 


Bo  1 


mo5^ 

The  dynamical  matrix  is  [fly]. 

Coupled  bending-torsion  vibration. — The  torsional  and 
bending  deflections  due  to  coupled  bending-torsion  vibra¬ 
tions  of  a  cantilever  beam  are  given  by  (appendix  B) : 


Ml 

To 


(12) 


where 


■  rgo^ 

1  lo  Bo 
5^  Co  Too 


ro  absolute  magnitude  of  projection  of  distance  from  elastic 
axis  to  center  of  gravity  on  perpendicular  to  bending 
direction  for  root  section 

r^o  radius  of  gyration  about  elastic  axis  at  root  section 
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The  quantities  and  are  defined  by  equations  (2) 
and  (8).  The  quantities  'yu  and  5<y  are  given  by 

S/cN'jii — {k — 1)  SicM' 'SrAf'jr”! 

k=l  L  r  =  k+l  J 

Sij  =  '^^-j^^(^Sk{iPjk—Qjk)-\-  +  y  (13) 

where 

{k—\y^jj{z)dz 

[j-l  (k-m+l  (k-iy'^fj{z)dz 

and  Sk  is  the  ratio  of  the  average  static  mass  unbalance  of 
the  A;***  interval  to  the  static  mass  unbalance  at  the  root 
section. 

The  Station  Numbers  P^k  and  are  listed  in  tables  I  to 
VIII  with  the  other  Station  Numbers.  The  determinantal 
equation  becomes 


Tail- 

X  r ci\2 

«ryu 

efiyia  . 

■  eUyiB 

r  a2i 

ra22  —  X  . 

•  k'a2„ 

erT2i 

€ry22  .  . 

•  eUyare 

ra„i 

ra„2 

•  •  r  Oinn~ 

-X  eUyni  ery„2  • 

•  'Ynn 

Sn 

5i2 

■  •  8in 

|Sn- 

X  ^12 

•  ^In 

821 

622 

■  •  82n 

^21 

^22  X  . 

■  ^2n 

8nl 

8n2 

.  .  ^nn 

fin2 

(14) 

or 

|X/-ho]|  =  0  (14a) 

where  [774^]  is  the  dynamical  matrix  and  I  is  the  identity 
matrix. 

The  first  n  roots  of  equation  (14)  will  give  the  first  n 
coupled  frequencies. 

APPLICATIONS  AND  RESULTS 

In  applying  the  previously  discussed  method,  it  is  necessary 
to  determine  for  a  given  beam  the  elements  oy,  /Sy,  yy,  and 
5y  of  the  dynamical  matrices.  These  quantities  wiU  depend 
on  the  physical  properties  of  the  beam  and  on  the  number  of 
stations  chosen.  If  the  physical  properties  of  the  beam  are 
known,  the  quantities  ay,  /Sy,  yy,  and  5y  can  be  directly  cal¬ 
culated  from  equations  (2),  (8),  and  (13).  The  numbers 
Mjk,  Njk,  Pjk,  Qjk,  M'ik,  N'jk,  P'jk,  and  Q' jk  appearing  in  these 
equations  depend  on  the  number  of  stations  n  that  are  used 
and  can  be  read  directly  from  tables  I  to  VIII  for  .any  given 
number  of  stations  up  to  eight.  Once  these  quantities  have 
been  calculated,  equations  (6),  (11),  or  (14)  can  be  solved 
for  the  frequencies  by  any  method  desired.  The  matrix- 


iteration  method  used  herein  is  simple  and  rapid  and  re¬ 
quires  no  particular  computing  skill.  As  will  be  indicated, 
however,  the  accuracy  of  equations  (6),  (11),  and  (14)  is 
such  that  relatively  few  stations  need  be  used,  in  which 
case  it  may  be  convenient  to  expand  the  determinants  and 
to  solve  the  resultant  low-order  algebraic  equation. 

In  order  to  illustrate  the  accuracy,  this  method  was  applied 
to  torsional  vibrations,  bending  vibrations,  and  coupled  vi¬ 
brations  of  a  uniform  cantilever  beam.  The  exact  theoretical 
values  for  torsional  vibrations  and  bending  vibrations  of 
uniform  cantilevers  are  well  known.  The  exact  theoretical 
values  for  the  coupled  bending-torsion  vibration  of  a  uniform 
beam  were  calculated  (appendix  D).  A  comparison  was 
then  made  between  the  values  obtained  by  the  method 
presented  and  the  exact  theoretical  values.  The  number  of 
stations  used  was  1,  2,  and  3  (n=],  n=2,  and  n=3).  The 
comparisons  are  summarized  in  table  IX. 

Torsional  vibration. — For  the  case  of  a  uniform  beam, 
Ck=Ik=l  and  equation  (2)  becomes 

Njk—{k  —  l)Mjk-\- (15) 

k  =  l  L  r=fc+l  J 

The  values  of  Njk  and  My*  are  given  in  tables  I  to  VIII. 
The  table  to  be  used  depends  on  the  choice  of  the  number 
of  stations. 


The  exact  theoretical  value  for  the  first  torsional  frequency 
is 

CO-1.571 

The  percentage  error  is  —1.4  when  only  one  station  is  used. 

The  mode  shape  obtained  by  the  method  of  Station  Func¬ 
tions  agrees  well  with  the  theoretical  mode  shape,  as  is  shown 
in  figure  2  (a). 

Let  71=2;  then  by  equation  (15)  and  table  II, 


«ii= A^ii+  Afi2 — Y5'^T2  — 


60 


AT  ,  Ar  31  ,  29  57 

ai2  —  N2i+  M22  240  "*"48  120 

AT  I  AT  8.8  16 

a2i — -Nil + — Jg -f  Yg  — 


AT  I  AT  31  ,  239 

0(22  — A/21 -I-M22—  240  ”^240 
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Pi  stance  from  root  x 
Length  of  beam  ’  X 


^d)  First  mode,  =3. 
(c)  Second  mode,  n=3. 
(0  Third  mode,  «=3. 


(a)  First  mode,  »=1. 

(b)  First  mode, 

(c)  Second  mode,  71=2. 

Figure  2.— Comparison  of  theoretical  mode  shapes  with  mode  shapes  obtained  by  taking  n  stations  along  the  beam  for  torsional  vibrations. 


The  determinantal  equation  then  becomes 
60  ^  120 

=  0 

16  13  ^ 

15  15“^ 

which  gives 

Xi=1.6214 

X2=0.1953 

Therefore 

"■=‘-W§ 

The  exact  theoretical  values  are 


0)1=1.571 

The  precentage  errors  of  the  first  two  modes,  for  only  two 
stations,  are  found  to  be  0  and  — 4.  ' 

The  mode  shapes  are  shown  in  figures  2  (h)  and  2  (c). 
Agreement  of  the  first  inode  with  the  exact  theoretical  shape 
is  excellent;  the  second  mode  agrees  fairly  well. 


Let  71=3;  then  by  equation  (15)  and  table  III, 

a„=iV„q-Mi3+Afi3=0.945833 
0112 = iV2i  “f"  Af22  “h  A/23  0 . 9 583 3 3 
a,3=iV3i+M32+A/33=0.520834 

a2i=iVu4-A^i2+2A/i3=  1.033333 
0)22”  A21“h  A^22“1~2A/23=  1 .883333 
0:23— ATsid"  A32“t"2A/33=  1 .011113 
a3i=Arii+Arj2+JVi3=  1.012500 
0)32 = ^21  “h  A^22  “h  ^23 = 2 .025000 
a33  =  iV31-|-iV32-l-Ai^38=  1 .387501 

The  determinantal  equation  is 

0.945833-X  0.958333  0.520834 

1.033333  1.883333-X  1.011113 

1.012500  2.025000  1.387501-X 

The  solutions  are 

Xi=3.6474 

X2=0.4093 

X3=0.1599 

Therefore 

«=4.689.V/#- 

„.-7.502.^^, 
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(a)  First  mode,  n=l.  (d)  First  mode,  w=s3. 

(b)  First  mode,  n=2.  (e)  Second  mode,  n=3. 

(c)  Second  mode,  n=2.  (0  Third  mode,  w=3. 


Figure  3.— Comparison  of  theoretical  mode  shapes  with  mode  shapes  obtained  by  taking  n  stations  along  the  beam  for  bending  vibrations. 


The  exact  theoretical  values  are 

.«=7.854^A 

The  percentage  errors  of  the  first  three  modes,  calculated 
by  use  of  three  stations,  are  found  to  be  0,  —0.5,  and  —4.5, 
respectively. 

The  mode  shapes  are  shown  in  figures  2  (d)  to  2  (f).  The 
first  two  modes  agree  very  well  with  the  theoretical  shapes ; 
agreement  of  the  third  mode  is  fair. 

This  procedure  can  be  carried  out  as  shown  for  any  number 
of  stations  desired. 

Bending  vibrations. — For  a  uniform  beam,  B^=mk=  1  and 
equation  (8)  becomes 

k=l(  r=fc+l  L\ 

(P_(^_1)3-(2^^  '  (16) 

Let  1 ; 

•  Pn—P'n —  Q'n 


and  from  table  I 

71  31  59 

630  1008  720 

Therefore,  from  equation  (7), 
co=3.493 

The  exact  theoretical  value  is 
w=3.516 

The  precentage  error  for  just  one  station  is  found  to  be 
-0.65. 

The  mode  shape  is  shown  in  figure  3  (a)  and  is  seen  to 
agree  very  well  with  the  theoretically  exact  shape. 

Let  n=2;  then  by  equation  (16)  and  table  II, 

^ii=P'ii-Q'ii+|iV'i2-^M'i2=0.422745 

0i3=F'2i-Q'2i+^iV'22-|M'22=O.295925 

^2i=2P'ii+2P'i2 —  Q'n —  Q' \2— 1.145 167 
^22  =  2P^21  4"  2P^22  —  Q' 2\  —  Q' 22=^  0.905530 
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The  characteristic  equation  is 


0.422745-X 

1.145167 


0.295925 

0.905530-X 


=  0 


The  roots  arc 


X,=  1.2943 


X2==0.0339 

The  exact  theoretical  values  are 

The  percentage  errors  for  two  stations  are  therefore  found 
to  be  0  for  the  first  mode  and  —1.5  for  the  second  mode. 
The  mode  shapes  are  plotted  in  figures  3  (b)  and  3  (c).  The 
first  mode  agrees  excellently  with  the  theoretically  exact 
shape;  the  second  mode  agrees  fairly  well. 

Let  n=3;  then  by  equation  (16)  and  table  III, 


Pn^P'u-Q'n+\N',,+^N\,-l  M\o-\  M'a3-0.270604 


The  characteristic  equation  is 
0.270604  — X  1.009943 
0.648170  3.266250  — X 

0.985135  5.822852 


0.487441 

1.689891 

3.204301-X 


The  roots  are 

X, =  6. 5521 
X2=0.1667 
X3=0.0223 

Therefore  _ 

a;i  =  3.516 

co2=22.04  ^^ 


The  exact  values  are 

0,1=3.516 

'0,2  =  22.04 


/3i2 — P'n —  A^^22”I“2 "^^23  —  g  M'22 — ^  g  M'ii —  1.009943 

/3,3  =  P'3,-e'3.  +  iA^'32  +  ^A^'33-^  M' M'33  =  0 .4  8  74  4 1 

021  =  27’'n-t-  2P'i2 —  Q' ii~  Q' \i-\~ 

|A^',2+2iV'i3-|M',2-|M'i3=0. 648170 

022=  2P'21-|-  2P'22—  Q'21—  Q^22  + 

|iV'22+  2iV'23-|  M'22-|  M'23=3. 266250 

023=2P^31“t"  2P'32 —  <2^31 —  Q'z2~{' 

I  fV'32  +  2iV'33  - 1  M'32  - 1  M'33  =  1 . 6  8  9  8  9 1 

031  =  3P'ii4  3P'i2+  3P'i3 —  Q'li —  Q'i2 —  Q! \i~\~ 

|Ar'i2+4A^'i3-^M',2-^M'i3=0.985135 

032=3P'21“["3P'22“I“3P'23 —  Q  21 —  Q' 22  —  Q'23'h 

I  A^'22  +  4Ar'23-^M'22-^M'23=5. 822852 

033=3P'3I  +  3P'32-1-3P'33 — Q' z\ — Q' i2 — ^'33+  ' 

|iV'32  +  4iV'33-|M'32— y  M'33  =  3.204301 


The  percentage  errors  for  three  stations  are  found  to  be  O, 
0,  and  —2.4,  respectively.  The  modes  are  plotted  in  figures 
3  (d)  to  3  (f).  The  first  two  modes  are  seen  to  agree  very 
well  with  the  theoretical  mode  shape;  agreement  of  the 
third  mode  is  fair. 

Coupled  bending-torsion  vibrations. — A  uniform  beam 
with  the  following  constants  was  chosen: 

7=^=38.56 

Ws 

6=0.8 


193.2 


er= 


The  values  of  an  and  are  obtained  as  previously  and  are 
the  same  as  given  before  for  n  =  l,  n=2,  and  n=5.  Also, 
because  iS*=Pj:=C)t=m4=4=l,  equations  (13)  become 


7ii=sr^4.-(fc-i)M',.+  i:  M'J 

k  =  l  L  r  =  k+l 


928716—51 - 2 
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Let  71=1;  then  the  determinant  is 


ra,i~X 

0.002156-X 

0.001196 

5ii 

^ii-X 

0.111111 

0.081944-X 

The  roots  are 

Xi=0.0837 

X2=0.0005 

"■==>•« 

The  procedure  for  calculating  the  exact  theoretical  values  is 
derived  in  appendix  D.  The  exact  values  are 

V5> 


The  percentage  error  for  the  first  mode,  calculated  by  use  of 
one  station,  is  —0.9. 


Let  71=2;  then  the  determinant  is 


ran  — X  raj2 


er7„  tPyja 


Taai 

6n 


ra22— X  «r72i  €r722 

^12  /3n — X  /3i2 


=  0 


s 


21 


^21  ^22 — X 


Substituting  the  known  values  and  solving  for  X  give  for  the 
first  two  roots 

Xi=1.3197 


X2=0.0412 

and  the  frequencies  become 

The  percentage  errors  for  two  stations  are  —0.3  for  the  first 
mode  and  —4.4  for  the  second  mode. 


This  procedure  can  be  carried  out  for  any  number  of  sta¬ 
tions  desired.  For  three  stations,  the  frequencies  obtained 

are 

/  £o 

3.48 

mol* 

co2=20.6 

Mol* 

co3=  48.2 

mjj* 

The  precentage  errors  are  —0.3  for  the  first  mode,  0  for  the 
second  mode,  and  —1.8  for  the  third  mode. 

The  results  obtained  by  the  method  presented  are  seen 
to  agree  very  well  with  the  exact  theoretical  values. 

These  results  are  summarized  in  table  IX,  where  a  com¬ 
parison  is  made  with  the  results  obtained  for  uncoupled 
bending  and  torsional  vibrations  by  use  of  influence  coeffi¬ 
cients  with  weighted  matrices  (reference  12).  The  values 
using  weighted  matrices  were  taken  from  table  I  of  refer¬ 
ence  12.  It  can  be  seen  that  for  a  given  number  of  stations, 
the  results  obtained  by  the  method  presented  herein  are  con¬ 
siderably  better  than  those  obtained  by  using  influence  co- 


Figure  4.— Variation  of  frequency  ratio  fi  with  coupling  coefficient  e  for  several  values  of 
uncoupled  frequency  ratio  T. 

efficients  with  weighted  matrices.  In  general,  it  is  indicated 
that  for  a  uniform  cantilever  beam  using  n  stations  along 
the  beam,  the  first  n—\  frequencies  and  modes  are  in  ex¬ 
cellent  agreement  with  exact  theoretical  values  and  even 
the  71®  mode  is  given  within  the  accuracy  with  which  the 
physical  properties  of  the  material  are  known.  For  a  tapered 
beam,  more  stations  may  be  required,  depending  on  the 
amount  of  taper.  The  number  of  stations  required  to  give 
satisfactory  accuracy  is  listed  in  table  X.  A  comparison  is 
made  by  using  weighted  influence  coefficients;  the  values 
are  taken  from  table  II  of  reference  12. 

The  first  vibrational  frequency  is  given  approximately  by 
equation  (C2)  (appendix  C)  when  coupling  exists  between 
bending  and  torsion;  it  is  plotted  in  figure  4.  In  order  to 
check  these  curves,  the  exact  solution  was  obtained  (appen- 
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dix  D)  for  the  ratio  {uitlostY  equal  to  4  and  Was  plotted  on 
the  same  figure.  The  values  given  by  equation  (C2)  are 
seen  to  be  in  excellent  agreement  with  the  theoretically 
exact  values. 

The  effect  of  the  coupling  between  bending  and  torsion 
is  to  reduce  the  first  natural  frequency  below  that  which 
would  exist  if  there  were  no  coupling.  This  effect  is  shown 
in  figure  4,  wherein  the  value  of  0  is  always  less  than  1. 
This  decrease  in  the  first  natural  frequency  due  to  coupling 
is,  however,  relatively  unimportant  in  the  practical  range  of 
(co,/co6)^>4  and  e<0.75. 

SUMMARY  OF  RESULTS 

A  method  based  on  the  use  of  Station  Functions  is  pre¬ 
sented  for  calculating  uncoupled  and  coupled  bending-torsion 
modes  and  frequencies  of  arbitrary  continuous  cantilever 
beams.  The  results  of  calculations  made  by  this  method 


indicated  that  by  the  use  of  Station  Functions  derived  herein, 
n  modes  and  frequencies  can  be  obtained  with  sufficient  ac¬ 
curacy  by  using  just  n  stations  along  the  beam  if  the  beam  is 
uniform.  For  a  tapered  beam,  more  stations  may  be  re¬ 
quired,  depending  on  the  amount  of  taper.  The  amount  of 
computational  labor  is  markedly  less  than  for  other  methods. 
The  use  of  Station  Numbers  tabulated  herein  further  re¬ 
duces  the  amount  of  calculation  necessary.  The  effect  of 
coupling  between  bending  and  torsion  is  shown  to  reduce  the 
first  natural  frequency  to  a  value  below  that  which  it  would 
have  if  there  were  no  coupling. 


Lewis  Flight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  October  18,  1949. 


APPENDIX  A 

SYMBOLS 


The  following  symbols  are  used  in  this  report: 

dij, 

coefficient  in  equation  for  Station  Function 
in  torsion 

B 

bending  stiffness  of  beam,  function  of  2 

Bo 

bending  stiffness  at  root  section  of  beam 

B, 

ratio  of  average  bending  stiffness  of  P*' 
interval  to  bending  stiffness  of  root 
section 

b  ij 

coefficient  in  equation  for  Station  Function 
in  bending 

C 

torsional  stiffness  of  beam,  function  of  2 

Co 

torsional  stiffness  of  root  section  of  beam 

c. 

ratio  of  average  torsional  stiffness  of  A:**' 
interval  to  torsional  stiffness  at  root 
section 

Cl,  C2,  C3 

constants  defined  in  appendix  B 

f}{z) 

Station  Function  in  torsion  for  j*"  station 
(defined  in  text) 

9}{z) 

Station  Function  in  bending  for  j*"*  station 
(defined  in  text) 

I 

mass  moment  of  inertia  per  unit  length  of 
beam  about  elastic  axis,  function  of  2, 
except  where  otherwise  defined 

h 

mass  moment  of  inertia  per  unit  length  of 
beam  about  elastic  axis  at  root  section 

h 

ratio  of  average  mass  moment  of  inertia 
per  unit  length  of  interval  to  mass 

moment  of  inertia  per  unit  length  at  root 
section 

i,j,k,n 

station  indices 

j,k,r 

summation  indices 

1 

length  of  beam 

Mjic,  Njic,  P jic, 

Station  Numbers  (defined  in  text) ;  function 

Qjk,  Ik, 

P'jk,  Q' Ik 

jk,  of  indices  j,  k,  and  n 

m 

mass  per  unit  length  of  beam,  function  of  2 

mo 

mass  per  unit  length  of  beam  at  root  section 

mk 

ratio  of  average  mass  per  unit  length  of 
A:**'  interval  to  mass  per  unit  length  at 
root  section 

n 

number  of  stations  along  beam 
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2j.(2) 

3<(2) 


rgo 


S 

So 

S, 


X 

y 

Vi 

z 

“tj!  ytji 
Vij 

r 

y 

5 


d 

di 

X 

n 

w 

W6 


bending  loading  function  on  beam 
torsional  loading  function  on  beam 
absolute  magnitude  of  projection  of  distance 
from  elastic  axis  to  center  of  gravity  on 
perpendicular  to  bending  direction 
radius  of  gyration  about  elastic  axis  at 
root  section 

absolute  magnitude  of  projection  of  distance 
from  elastic  axis  to  center  of  gravity  on 
perpendicular  to  bending  direction  for 
root  section 

static  mass  unbalance,  function  of  z,  mr 
static  mass  unbalance  at  root  section,  moro 
ratio  of  average  of  static  mass  unbalance  at 
it**'  section  to  static  mass  unbalance 
at  root  section 

distance  from  root  of  beam,  except  where 
otherwise  defined 
bending  deflection,  function  of  z 
bending  deflection  at  1“*  station 
dimensionless  distance  along  beam,  xjb 
elements  of  dynamical  matrix  defined  in  text 

5^  Co  mo 

uncoupled  frequency  ratio, 
length  of  interval  along  heam  between 
two  stations 

coupling  coefficient,  (roIrgoY 
torsional  deflection,  function  of  2 
torsional  deflection  at  i**'  station 
root  of  frequency  equation  or  characteristic 
root  of  dynamical  matrix 
frequency  ratio,  (w/co^)^ 
frequency  of  vibration 

frequency  of  uncoupled  fundamental  bend¬ 
ing  mode 

frequency  of  uncoupled  fundamental  tor¬ 
sional  mode 

second  derivative  of  deflection  with  respect 
to  time 


m 


APPENDIX  B 


STATION  FUNCTIONS  AND  DETERMINANTAL  EQUATIONS 


TORSIONAL  VIBRATIONS 

A  schematic  diagram  of  a  cantilever  beam  divided  into  n 
intervals  of  length  5  is  shown  in  figure  1.  The  Station 
Functions  for  the  torsional  vibrations  of  such  a  beam  must 
satisfy  the  following  conditions: 


2  =  0 

/do)=o 

(Bl) 

z=n 

ft(n)  =  0 

(B2) 

z=i 

/i(i)  =  l 

(B3) 

z=j 

fi(j)  =  0  j9^i 

(B4) 

where /'(2)  denotes  the  derivative  with  respect  to  z. 

Equations  (Bl)  and  (B2)  represent  the  boundary  condi¬ 
tions  that  must  be  satisfied  by  a  cantilever  beam  vibrating 
in  torsion;  equations  (B3)  and  (B4)  represent  the  further 
conditions  imposed  upon  the  Station  Functions.  These 
conditions  will  be  satisfied  by  a  function  of  the  type 

/i(2)=au2+a2,2:^+  .  •  •  +a(«+i)i2'’‘+‘'  (B5) 

where  the  coefficients  a,,  must  satisfy  the  following  simul¬ 
taneous  equations  obtained  from  conditions  (B2),  (B3), 
and  (B4): 


0  —  .  .  .  -(-(nT l)u-”<i(K+i)i  (B2a) 

\  =  .  .  .  -|-f  ci  (B3a) 


0=iaii-fi^a2i4-i®a3i+  .  .  .  (B4a) 


The  coefficients  can  be  obtained  by  solving  equations 
(B2a)  to  (B4a)  and  the  functions  jt{z)  determined  for  each 
station.  Equation  (B5),  however,  can  also  be  written  in 
the  following  form: 


/i(z) 


n(2— i)2(2-Ci) 

_ 

^  — Cl) 

jVi 


(B5a) 


where  11  represents  the  product  for  all  values  of  j  except 

j=i.  The  function  in  equation  (B5a)  obviously  satisfies 
conditions  (Bl),  (B3),  and  (B4)  because  it  has  zeros  at  all 
points  specified  by  equation  (B4),  it  equals  1  at  the  point 
specified  by  equation  (B3),  and  it  equals  zero  at  the  point 
specified  by  equation  (Bl).  In  order  to  satisfy  condition 
(B2),  the  constant  Ci  is  determined  by  substitution  of  equa¬ 
tion  (B5a)  into  equation  (B2). 

Ci=n  for  i 


Ci=n 


( 


1+- 


1+2: 

jj^n 


for  i  =  n 


Equation  (B5)  can  be  obtained  from  equation  (B5a)  by 
carrying  out  the  indicated  multiplications.  The  complete 
deflection  function  is  then  given  by 

®(2)=/l(2)®l+/2(z)6'2+  ■  •  ■  +/«(2)0« 

=  (B6) 

The  continuous  loading  function  g*  (s)  can  now  be  written 
as 

qt{x)  =  Io}^6{z)  =  Iw'^'^jj{z)6j  (B7) 

A  continuous  loading  function,  which  is  a  function  of  the 
deflections  at  the  reference  stations,  has  thus  been  obtained. 

BENDING  VIBRATIONS 

The  Station  Functions  for  the  bending  vibrations  of  the 
beam  shown  in  figure  1  ftiust  satisfy  the  following  conditions : 
at 


2  =  0 

P4(0)  =  0 

(B8) 

2  =  0 

g'ii0)  =  0 

(B9) 

z—n 

g'\{n)^0 

(BIO) 

z=n 

g"\(n)  =  0 

(Bll) 

z  =  i 

fi'i(i)  =  l 

(B12) 

z=j 

<7i(i)  =  0  jVf 

(B13) 

where  g'  (2),  g"  (z),  and  g'"  (2)  denote  the  first,  second,  and 
third  derivatives,  respectively,  of  g  (2)  with  respect  to  2. 

Equations  (B8)  to  (Bll)  represent  the  boundary  condi¬ 
tions  that  must  be  satisfied  by  a  cantilever  beam  vibrating 
in  bending  and  equations  (B12)  and  (Bl3)  represent  the 
additional  conditions  imposed  upon  the  Station  Functions. 
These  conditions  will  be  satisfied  by  functions  of  the  type 

£^>(2)=^>2i2^+^>3i2^+ •  •  ■ +6u+3)i2'"'^^'  (B14) 

where  the  coefficients  bij  must  satisfy  the  following  equa¬ 
tions  obtained  from  conditions  (BlO)  to  (Bl3) : 

0  =  262i  +  6ii  63i+ .  .  . -{-("U-t- 3)('*^+2)U'^”"^** (B4-3)i  (BlOa) 

0  =  6f>3i -|- 2471644 -|- .  •  ■ (n+s)i  (Blla) 

1  =  7^62i++634-h  .  .  • 6 („+3)i  (Bl2a) 

0=+62i-i-+63i+ .  .  . -f  i'''+®’6(„+3)f  jy^i  (Bl3a) 

The  coefficients  can  therefore  be  obtained  from  equations 
(BlOa)  to  (Bl3a)  and  the  fimctions  gt  (z)  determined  for 
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each  station  i.  Equation  (B14)  can,  however,  be  written 
in  the  following  form: 

U{z—j)z^{z^  +  C2Z-\-Ci) 

where  n  represents  the  product  for  all  values  of  j  except 

jVi 

j—i.  The  function  in  equation  (Bl4a)  obviously  satisfies 
conditions  (B8),  (B9),  (Bl2),  and  (B13),  because  it  has 
zeros  at  all  points  specified  by  conditions  (B8),  (B9),  and 
(B13)  and  equals  1  at  the  point  specified  by  equation  (B12). 
In  order  to  satisfy  conditions  (BlO)  and  (Bll),  the  constants 
Ca  and  C3  are  determined  by  substitution  of  equation  (Bl4a) 
into  equations  (BlO)  and  (Bll).  The  general  forms  for  Ca 
and  C3  are,  however,  complicated  and  it  is  easier  to  obtain 
the  numerical  values  of  these  constants  for  each  specific 
case.  Equation  (B14)  can  then  be  obtained  from  equation 
(Bl4a)  by  carrying  out  the  indicated  multiplications.  The 
complete  deflection  function  is  then  given  by 

yis)='Ilgiis)yi  (B15) 

j=t 

The  continuous  bending  loading  function  q^^z)  can  now  be 
written  as 

qtiz)^mw^y{z)=-moi^J29ii2)yi  (BI6) 

COUPLED  BENDING-TORSION  VIBRATIONS 

The  Station  Functions  for  the  coupled  bending-torsion 
vibrations  are  the  same  as  previously  given  for  the  bending 
vibrations  and  the  torsion  vibrations.  The  loading  func¬ 
tions,  however,  are  given  as  follows  (reference  7): 

qi{z)  =  Ioi^d(,z)  +  S(Jy{z) 

=  <^^'i2lIJjiz)0j+Sgj{z)yj]  (Bl7) 

;=i 

and 

q„{z)=S  00^0(2)  -f  y{z) 

=  0)^2  [Sfjiz) di + mgjiz)  yj]  (B 1 8) 

DETEHMINANTAL  EQUATIONS  AND  DYNAMICAL  MATRICES 

Once  the  Station  Functions  and  the  corresponding  loading 
fimctions  have  been  determined,  the  deflections  at  the 
reference  stations  can  be  obtained  in  terms  of  the  loading 
fimction.  A  homogeneous  equation  in  the  reference-station 
deflections  for  each  station  is  thereby  obtained.  The 

determinant  of  the  coefficients  of  the  resultant  set  of  homo¬ 
geneous  equations  can  be  set  equal  to  zero ;  the  determinantal 
frequency  equation  is  thus  derived.  The  deflections  at  the 
reference  stations  are  obtained  by  the  well-known  equations 
for  obtaining  influence  coefficients. 

Torsion. — The  deflection  at  the  station  i  due  to  the 
continuous  loading  qt{z)  on  the  beam  is  given,  by? ' 


=  ^  dz  (B19) 

If  C  is  assumed  to  have  a  constant  value  for  each  interval, 
these  integrals  may  be  written  as  the  sum  of  integrals  over 
each  section.  Equation  (B19)  then  becomes 

zqt{z)dz-\-{  {l—k)q,{z)dz+f  g, (2)^21 
Oo  fc=l  —  l  J  k—1  Jk  J 

(B20) 

By  substituting  the  relation 

qt{z)  =  u)U  Z;/;(2)0j 
;=i 

and  by  assuming  a  constant  value  for  I  for  each  interval  and 
changing  the  summation  order, 

6»i = 0)^52^  S  I S  7^  f  zfiiz)  dz—{k  —  l)lA  Jj{z)dz + 
Uo;  =  l  f*=lU*:L  J  k-l  Jk-1 

(B21) 

Let 

rk 

zji{z)dz=Nj„ 

ct'  \  (B22) 

fj{z)dz=Mj^ 

Jk-l  J 


where 


9i=  Op’  -7^  P  Xj  Olil 
Uo  j=l 


y=r  [~ Ij  Nile. — {k  1)  I/cMjic-]-  (B24) 

k=l  L  r=t+l  J 


If  Ck=Ik=l  (constant  cross  section),  then 


Uij — ^Jk — (k  —  A/jr”|  (B25) 

*=1  L  r=t+l  J 


;=i 

and  the  characteristic  equation  is 


(B23a) 


l[a,^]-X7|  =  0  (B27) 

where  I  is  the  identity  matrix. 

Bending. — The  deflection  at  the  station  i  due  to  the  con¬ 
tinuous  loading  qt,(z)  on  the  beam  will  be  given  by 


gi.(3)f 

Jo  Jo 


(2  — 2i)(-i  — 2i) 


/•«  pi  r?  — 


(2  — 2i)(f  — 2i) 


dzidz-{- 


dzidz 


(B28) 
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If  B  is  assumed  to  have  a  constant  value  for  each  interval, 
these  integrals  may  he  written  as  the  sum  of  integrals  over 
each  interval.  Equation  (B28)  then  becomes 

jk-i  {k  —  lfz+~  (A:  — 1)=*J  (i^{z)dz  + 

ij^  |^2—|(2fc  —  l)Jg^6(2)cfs-l- 

(B29) 

By  substituting  the  relation 

qbiz)  =  io^m'^gj(z)yj  (B30) 


and  by  assuming  a  constant  average  value  for  m  in  each  in¬ 
terval  and  changing  the  summation  order, 


where 


n 


^hVi 

;=i 


(B31) 


ik~Q'  ik)-\-  S  mr[(i-k+^')N'jr  + 

k=\  Ok\  r  =  /i+l  L\ 


( 


{2k— \)  . 


)mv 


(B32) 


gjiz)  dz 


C'/t— ^  (^■— 1)^2+|  {t— l)d  g,{z)  di 

^'jk=  r  zg]{z)dz 
Jk-l 


(B33) 


M'}k=  r  gj{z)  dz 

Jk-l 

For  a  uniform  beam,  mk=Bk=l  and  equation  (B32)  be¬ 
comes 

fiii=iz(iP'ik-Q:ik+''iz  \(i~k^-^ 

k  =  l  \  r  =  k+l  (\  2, 

)iVV  + 

(B32a) 

Let 

\ _  Bo 

w^b^mo 

(B34) 

then  the  characteristic  equation  becomes 

T 

o 

(B35) 

where  7  is  the  identity  matrix  and  jSy  is  the  dynamical  ma¬ 
trix.  In  expanded  form,  equation  (B35)  becomes 


/3u-X 

012 

•  •  01n 

^21 

022  X  . 

■  02n 

0n2  •  ■ 

•  (^nn  ^ 

(B35a) 


where  X  is  a  latent  root  of  the  matrix  [j8ij]. 


Coupled  bending-torsion  vibrations. — The  deflections  at 
station  i  are  given  as  before  by  equations  (B19)  and  (B28). 
The  loading  functions  g,  and  gj  are  changed  as  follows: 


qi(z)  =  u^[1 0(z)  +  Sy(z)]) 
qt{z)  =  e{z)  -f-  m  y{z)\ ) 


(B36) 


If  these  two  equations  are  substituted  into  equations  (Bl9) 
and  (B28)  and  the  integrations  are  performed  as  previously, 
the  following  relation  is  obtained; 


(B37) 


where  ay  and  ,Sy  are  given  in  equations  (B24)  and  (B32)  and 

e=^ 

P _  1  7  c  j?o 

~  Coma 

S  ^  r  SkN'ik- {k-\)  SkM',k+  i:  Sr  M' J 
fk^-{k-\f  27-l.\,,  D 

where 

{k-lf'jjA^)dz 

the  deterniinantal  equation  therefore  is 
|X7— [r)y]|  =  0 

where  [r/y]  is  the  dynamical  matrix,  the  elements  of  which 
are  as  indicated  in  equation  (B37).  The  matrix  [rjy]  is  seen  to 
be  a  2nX  2n  matrix. 


APPENDIX  C 

QUADRATIC  FORMULA  FOR  FIRST  COUPLED  MODE 


If  only  the  first  vibrational  mode  is  desired,  it  is  possible 
to  obtain  this  mode  approximately  by  coupling  together  the 
fundamental  uncoupled  bending  mode  with  the  fundamental 
uncoupled  torsional  mode  to  obtain  a  simple  quadratic 
equation  for  the  first  coupled  frequency.  This  equation  is 
valid  when  the  coupling  coefficient  e  is  constant  along  the 
beam.  The  differential  equations  obtained  by  coupling  the 
fundamental  uncoupled  torsional  mode  with  the  funda¬ 
mental  uncoupled  bending  mode  are: 

my + 6  +  m  = 0 

8y 16 1 u3t^d=Q  3 

where 

m  mass  per  unit  length  of  beam,  function  of  z 
S  static  mass  unbalance,  function  of  z 


I  mass  moment  of  inertia  about  elastic  axis,  function  of  2 
cos  frequency  of  uncoupled  fundamental  bending  mode 
o>t  frequency  of  uncoupled  fundamental  torsional  mode 
. .  denotes  differentiation  twice  with  respect  to  time 

These  equations  lead  to  a  quadratic  equation  in  the  fre¬ 
quency  ratio  11,  whose  solution  for  the  lowest  frequency, 
provided  e  is  constant  along  the  beam,  is 


where 


1—7 

2(1 -t) 


47(1  —  «)' 

(1_,A2 


(C2) 


H  frequency  ratio, 

7  uncoupled  frequency  ratio, 

€  coupling  coefficient,  (r/rgY 

This  quadratic  has  been  plotted  in  figure  4  for  values  of  e 
ranging  from  0  to  1  and  values  of  7  =  (utloitY  from  1  to  100. 


APPENDIX  D 

EXACT  SOLUTION  FOR  COUPLED  BENDING-TORSION  VIBRATIONS  OF  UNIFORM  CANTILEVER  BEAM 

where 


The  differential  equations  for  the  equilibrium  of  an  ele¬ 
ment  of  a  beam  vibrating  in  coupled  bending-torsion  vibra¬ 
tions  can  be  put  in  the  following  dimensionless  form: 


d*Yi  ml^  .  ml^ 

ijF-Tr"^'+Tr"'‘ 

d^Y2 


IV  IF 

dx^  *  (7  “  C 


where 


Now 


Yi^ylr 

Y2=e 

distance  from  root 

'  I 

t^ir/rgY 

2  CiB 


(>>t  —Cs 


IF 


where 


C4=12.36 

C5=2.467 

Equations  (Dl)  become 

d^Y, 


dx^ 

d^_ 

dx^ 


CiQ{Y  id-  T2) 


:-6  —  F2 

7  7 


(Dl) 


(D2) 


Let 


Then 


(fF, 

dx 

dY, 

dx 

dY, 

dx 

dY, 

dx 

dY_, 

dx 


^6 

dx 


il^iwIwiY 

7=(oj,/co6)^ 

=  F3 

=  Y, 

=  F, 

=  Fe 

=C4ll(Fi  +  F2) 


(eF,+  F2) 


7 


(D3) 


Equation  (D3)  can  be  written  as  the  single  matrix  equation 


A 

dx 


“Fr 

“  0 

0 

1 

0 

0 

0“ 

“Fi" 

F2 

0 

0 

0 

0 

0 

1 

F2 

F3 

0 

0 

0 

1 

0 

0 

F3 

F4 

= 

0 

0 

0 

0 

1 

0 

F4 

F, 

Ci9, 

C4D 

0 

0 

0 

0 

Fa 

Fa 

—  €  C5O 

0512 

0 

0 

0 

0 

Fe 

_  7 

7 

(D4) 
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IS 


ax 


(D4a) 


where  Y  and  A  are  the  matrices  indicated. 

The  solution  to  the  matrix  equation  (D4)  is  given  by 

Y=e^^Y, 

where  Fo  is  a  column  of  arbitrary  constants. 

From  the  boundary  conditions 


(D5) 


at  x  =  Q 
x=l 


Fi=F2=F3=0 

F4=F3=Fe=0 


Fo=F(0)  = 


0 
0 
0 

F4(0) 
LY,{0)J 


If  then  Qij  is  an  element  of  the  matrizant  e^,  the  boundary 
conditions  give 


O44 

064 


^45 

^^55 


046 

^56 

Osfi 


=  0 


(D6) 


Equation  (D6)  is  the  frequency  equation.  It  has  an  infinite 
number  of  roots  for  co. 

In  order  to  determine  the  elements  Oj^,  must  be  evalu¬ 
ated.  Use  will  be  made  of  Sylvester’s  theorem  (reference 
13). 


The  X  matrix  of  A  is 


0 

0 

0 

C4O 

C5O 


0 

—  X 
0 
0 

C40 

C50 

7 


0 

0 

1 

-X 

0 

0 


0  ■ 
1 
0 
0 
0 

-X 


L  7 

The  characteristic  equation  A(X)  =  0  is 

CsO 


02 


X«+^  X"-C40X2-(l-e)C4C5  — =0 
7  7 


(D7) 


Equation  (D7)  is  a  cubic  equation  in  X^.  Let  the  roots  be 

^2,— ^2)  ^3,-^3 

Then  by  the  confluent  form  of  Sylvester’s  theorem, 

=  _ ^  r — — 1  mg) 

^|(ai-l)dX“-i  n(X-X,)“^  ,  ^ 

' -  Lkr^i  = 

where  F(\)  is  the  adjoint  matrix,  r  is  the  number  of  distinct 
roots,  and  a*  is  the  multiplicity  of  the  F*  root. 

If  the  roots  are  all  distinct,  this  relation  becomes 


e 


e\F{\t)—e  \F(— Xj) 
fif  2X4n(X4-X^)(X4  +  X^) 


(D9) 


where  the  adjoint  matrix  E(X)  is  given  by 


X«+^  X3 

7 

C40X 

X^-h^  X2 

7 

x^-t-^x 

7 

X^-hCs- 

7 

C4O 

-tCs-X' 

7 

X®— C4OX 

-6^  V 

7 

C5O 

—  e  —  X 

7 

Csfi 
—  e  — 

7 

X^— C4O 

02 

C4fiX^+  (1  —  e)c4C5  — 

7 

C40X'* 

X=+^  x^ 

7 

x^-t-^  X^ 

7 

x^-h^x 

7 

(/4t2X 

Q2 

C4fiX^+(l  —  €)C4C5  -  X 

7 

C40X^ 

C40X2+(1  e) 

7 

X«  +  ^  X' 

7 

X4-(-^  X' 

7 

C4fiXH(l— e)C4C6  — X=“ 
7 

C4nx* 

C4i2XH(l  e)‘^‘^*^'x 

7 

C4JlX^+(l-f)  — 

7 

X*+^  x^ 

7 

c,Q\^ 

G\ 

1 

1 

_ 1 

7  7 

-6^X^ 

7 

7 

CsQ  . 

—  e  —  X 

7 

X'— C4fiX 

(DIO) 
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From  equations  (D9)  and  (DlO),  the  elements  are  seen 
to  be  given  by 

\  4  I  ^6^  N  2 

3  - At 

n  coshXt 


^4S —  — 


Xt'+'^  Xt 


tef  XtHCXt^-x/) 

i^i 


sinh  X< 


o  _  -A  CiUXi^  .  , 

^  Xt  n  (Xi^-  X/) 

j^i 

3  c40Xt^  +  ^(l-e)  _ 

“  S  Xtn(Xt'-x/) 

^^55=^44 

^  C4fiXt^ 

§  n(Xt^-x^^ 


5284= -S 


52^ 

—  eCs  —  A 


sinh  Xt 


n  (Xi^-x^^ 

52  ^ 

3  —  €C5  —  Xt 

Xtn(xAx,^ 

r.  ^  X/-C452 

§  n(Xt'-x,^ 

jVi 


(Dll) 


The  value  of  the  determinant  in  equation  (D6)  must  be  plot¬ 
ted  against  the  frequency;  the  value  of  the  frequency  for 
which  this  determinant  becomes  zero  is  thereby  obtained. 
This  procedure  involves  first  solving  the  cubic  equation  (D7) 


for  each  assumed  value  of  frequency  parameter  and  then 
calculating  the  elements  of  the  determinant  from  equations 
(Dll).  The  process  is  evidently  long  and  laborious. 
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TABLE  I— STATION  NUMBERS 


TABLE  II— STATION  NUMBERS 


A 

1 

2 

i 

M 

11 

5 

1 

12 

12 

N 

8 

8 

15 

15 

P 

0. 183333 

0. 025000 

Q 

M' 

.046032 

. 029365 

.536364 

. 627273 

N' 

.  367100 

.  851948 

P' 

.137933 

. 057955 

Q' 

.  036616 

.  069733 

M 

13 

29 

2 

48 

48 

N 

31 

239 

240 

240 

P 

-0. 037500 

0. 143750 

0 

-.008135 

.  181448 

M' 

-.060795 

. 448674 

N' 

-.034875 

. 758685 

P' 

-.011252 

. 118462 

Q' 

002614 

. 150415 
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TABLE  III— STATION  NUMBERS 


n=3 


ft 

i 

1 

2 

3 

M 

0. 950000 

0.  450000 

-0.  050000 

N 

,  546833 

.  587500 

-.120833 

P 

.  186310 

. 032143 

-.005357 

s, 

.  046577 

.  038244 

-.011756 

M' 

.  596268 

. 533205 

-. 097426 

N' 

.  399646 

.  708205 

-.239994 

P' 

. 148013 

.  042560 

-.012798 

Q' 

.  038884 

.  050843 

-.028318 

M 

2 

“0.  525000 

0.  725000 

0. 475000 

N 

241667 

1.175000 

1. 091667 

P 

068452 

.  160714 

. 031548 

Q 

014583 

.202083 

.  068750 

M' 

149356 

.602896 

.  625418 

N' 

-.083406 

.  994860 

1.  475163 

P' 

026378 

.  143948 

.  057937 

Q' 

006034 

.  181698 

.  127659 

M 

3 

0.  235185 

-0. 153704 

0.  568519 

N 

.  106019 

231944 

1.  613426 

P 

.029563 

-.023677 

.  139749 

fi, 

. 006222 

-.028963 

. 316408 

M' 

.  040630 

-.072928 

.  445812 

N' 

.  022325 

-.111744 

1.  200133 

P' 

.  006972 

012081 

.118007 

Q' 

.  001579 

-.  014830 

.  267865 

TABLE  V— STATION  NUMBERS 


w  =  5 


1 

2 

3 

4 

5 

M 

1 

1. 097991 

0. 408755 

-0. 040898 

0. 019866 

-0.013120 

N 

.  608222 

.  527493 

-.100112 

. 069159 

-.058445 

P 

.  202887 

.026908 

-. 005074 

.002776 

-.001627 

Q 

. 049943 

.031910 

-.011210 

. 008927 

-. 006836 

M' 

. 649902 

.  492141 

070298 

. 034939 

-.024007 

N' 

. 427616 

.  647530 

-.  172488 

.121519 

-.107639 

P' 

. 156411 

.  036908 

-.008903 

. 004824 

-.003376 

<?' 

.040729 

.  043977 

-.019678 

. 015509 

-.014228 

M 

2 

-0.  839550 

0.  799339 

0.  493783 

-0. 089560 

0. 049339 

N 

-.373049 

1.  282044 

1. 145470 

-.310549 

.  219544 

P 

-.103119 

.  169599 

.  037952 

-.011949 

.  006007 

Q 

-.021583 

.  212792 

.  083245 

038398 

. 025243 

M' 

-.255330 

.699256 

.523828 

-. 099723 

.  058134 

N’ 

-.  139170 

1. 138472 

1.219101 

345551 

.260447 

P' 

-.043239 

.157833 

. 041704 

-.013166 

. 008078 

O' 

009758 

.  198610 

. 091532 

042299 

. 034055 

M 

3 

0. 762798 

-0. 313591 

0. 651687 

0.  575298 

-0. 126091 

N 

.  329315 

-. 465823 

1. 718204 

1. 923065 

559573 

P 

. 089079 

-.044602 

. 150488 

. 049347 

-.014691 

0 

. 018a34 

-.054326 

. 340126 

.  157843 

-.061693 

M' 

. 197103 

-.228783 

.633812 

.  573549 

-. 137821 

N' 

. 105126 

-.344915 

1. 674943 

1. 916360 

-.616234 

P' 

. 032115 

-.034985 

. 148520 

.  048803 

-.018600 

Q' 

. 007151 

-.042759 

. 335798 

.  156076 

-.078385 

M 

4 

-0. 548214 

0. 187897 

-0. 159326 

0. 576786 

0.  562897 

N 

-.  233780 

.  276637 

-.400446 

2. 109970 

2. 4.32887 

P 

-.062665 

.  025479 

-.024868 

.140460 

. 042295 

0 

-.  012809 

.  030950 

-.055302 

.458397 

. 177006 

M' 

-.  117990 

. 117132 

-.  136462 

. 567359 

.664660 

N' 

062435 

. 175188 

-.344108 

2. 041363 

2. 901646 

P' 

-.018958 

.  017186 

021868 

. 137234 

.  063512 

Q' 

004201 

.  020954 

-.048664 

.447980 

.  267043 

M 

5 

0. 214238 

-0, 069026 

0. 050904 

-0. 080137 

0. 525349 

N 

.090928 

-.  101362 

.  127410 

-.283759 

2.  458336 

P 

.  024289 

-.009225 

.  007667 

-.013645 

.  134478 

0 

.  004952 

-.011196 

.  017031 

044065 

.  573757 

M' 

.  033722 

-.031711 

.  032307 

-.056459 

.  432107 

N' 

.  017789 

-.047311 

.081140 

-  200039 

2. 030260 

P' 

.  005389 

-.004593 

.  005002 

-.009675 

.  116059 

O' 

.001192 

-.005596 

.011120 

-.  031249 

.495663 

TABLE  IV— STATION  NUMBERS 


n~4 


■a  / 

/ 

1 

2 

3 

4 

M 

1 

1.022222 

0. 429630 

-0. 051852 

0. 022222 

N 

. 576455 

. 557937 

-.  127249 

. 076455 

P 

. 194478 

. 029597 

-.006581 

.002612 

0 

.  048240 

.  035167 

-.014647 

. 008359 

M' 

. 623188 

.511882 

-.082891 

. 042276 

N' 

. 413738 

. 676680 

-.203719 

. 146954 

P' 

. 152256 

. 039616 

-.010651 

.  005795 

Q' 

. 039818 

. 047267 

-.023551 

.018630 

M 

2 

-0. 647917 

0. 747917 

0. 618760 

-0. 085417 

N 

-.292857 

1. 207143 

1.  207143 

-.292857 

P 

-.081920 

. 163021 

. 041295 

-.  009598 

0 

-. 017295 

. 204828 

.090642 

-. 030688 

M' 

-.211987 

. 667412 

.544025 

112648 

N' 

116662 

1.091462 

1.269193 

-.390585 

P' 

-.  036502 

. 153469 

. 044508 

-.015000 

Q' 

-.008281 

. 193310 

. 097745 

-.  048203 

M 

3 

0.  522222 

-0.  255556 

0. 633333 

0.  522222 

N 

.  229365 

-.381746 

1.673810 

1.729365 

P 

.  062798 

037401 

.  148512 

.037202 

0 

.  013040 

045624 

. 335791 

. 118397 

M' 

. 122052 

166738 

. 582158 

. 643846 

N' 

. 065879 

-.  252823 

1. 545802 

2. 164827 

P' 

.  020304 

-.026235 

.  140822 

.060564 

Q' 

. 004551 

-.032114 

. 318707 

.  194016 

M 

4 

-0.  221701 

0. 094850 

-0. 105961 

0. 543924 

N 

-.096544 

. 140724 

-.267841 

1.997206 

P 

-.026267 

. 013391 

-.017322 

.136803 

0 

-.005428 

. 016301 

-.038574 

,  446729 

M' 

-. 035456 

. 042628 

-.064723 

. 438962 

N' 

-.019023 

. 064205 

164169 

1. 622066 

P' 

-.005836 

. 006481 

-.010869 

.117037 

Q' 

-.001303 

. 007917 

-.024222 

.  382752 

TABLE  VI— STATION  NUMBERS 
n=6 


j 

1 

2 

3 

4 

5 

6 

M 

1 

1. 172073 

0. 391101 

-0. 032371 

0.  013323 

-0. 010149 

0. 008879 

N 

.  638800 

.  501856 

-. 078978 

.  046300 

-.045644 

. 048522 

P 

.  210893 

.  024685 

003894 

.  001823 

-.  001498 

.001143 

Q 

,  051551 

. 029221 

008595 

.  005860 

-.  006322 

.  005949 

M' 

.  676394 

. 474177 

-.059129 

.  026682 

-.  018685 

.015649 

N' 

.  441269 

.  621067 

144759 

.  092370 

-.  083901 

.  085903 

P' 

. 160476 

.  034473 

-. 007337 

.  003631 

-.  002691 

.002241 

Q' 

.  041616 

.041021 

-.016206 

.  011674 

-.011350 

.011694 

M 

2 

-1.066598 

0.  853106 

0.  468124 

-0. 070505 

0.  044513 

-0. 035782 

N 

-.466718 

1.  360105 

1.  081893 

244062 

.  199948 

195451 

P 

127634 

. 176358 

.034411 

-.  009203 

.006452 

-. 004561 

Q 

-.  026505 

.  220969 

. 075401 

029565 

.  027216 

023740 

M' 

.303948 

.  731991 

. 503742 

085145 

.  049793 

-.  038661 

N' 

-. 164215 

1. 186681 

1. 169252 

-.  294736 

.  223322 

212149 

P' 

050692 

. 162263 

. 038896 

-.011103 

.  007043 

— .  005503 

Q' 

-.011384 

.  203987 

. 085309 

036665 

.  029701 

028708 

M 

3 

1. 150584 

-0.  404200 

0.  693794 

0.  546418 

-0. 133366 

0.  089627 

N 

.  489124 

597290 

1. 822457 

1.822457 

-.597296 

. 489124 

P 

. 130870 

-.055956 

. 156259 

.  045293 

-.  018484 

.011225 

Q 

.  026719 

-. 068060 

. 352909 

. 144808 

-.077925 

.  058410 

M' 

.  267118 

275585 

. 662165 

. 553477 

-.  126314 

.083246 

N' 

.  141177 

-.  413819 

1.  745286 

1.  846431 

-.  564890 

. 456509 

P' 

.  042839 

041308 

. 162473 

.  046980 

-.017100 

.011713 

Q' 

.  009490 

050434 

.  344556 

.  146998 

-.  072063 

. 061097 

M 

4 

-0.  930965 

0. 273028 

-0. 194854 

0.  592473 

0.  624591 

-0. 171416 

N 

390902 

. 399897 

488124 

2. 1^3786 

2.  720209 

-. 933437 

P 

103635 

.  036020 

-.  029594 

. 142229 

.  056698 

-. 020561 

Q 

-.021011 

. 043690 

-.  065761 

.  464052 

.  238215 

106938 

M' 

210538 

. 182630 

-. 180822 

. 598464 

.604424 

-.  166643 

N' 

-.110263 

.  271857 

-. 454522 

2. 185427 

2.  628772 

-.912352 

P' 

-.  033225 

.  026154 

-. 028221 

. 143453 

.  053359 

-.  022768 

Q' 

007320 

.  031846 

062753 

.468015 

.  224100 

-. 118729 

M 

5 

0. 581796 

-0. 156399 

0. 092907 

-0. 111954 

0.  537351 

0. 599157 

N 

.  242612 

228221 

.  231501 

394888 

2.  509279 

3. 194001 

P 

.063996 

-.  020218 

.  013474 

-. 018261 

.  134578 

. 046991 

Q 

,  012925 

024496 

.  029896 

058920 

.  674036 

.  243745 

M' 

.  120308 

-.097119 

.  081568 

110987 

.535339 

.  685021 

N' 

.  062735 

144101 

.204039 

-.  391731 

2.  499695 

3.  678582 

P' 

.018841 

013674 

.  012212 

018226 

.  133988 

.  066470 

Q’ 

.  004140 

-.016634 

.  027119 

-.  058815 

.  571521 

.  345986 

M 

6 

-0.  209220 

0.  054246 

-0. 030239 

0. 031561 

-0.  064035 

0.  510543 

N 

086982 

. 079042 

075223 

. 110987 

-.  291427 

2.  902746 

P 

-.  022893 

.  009657 

004323 

.004969 

-.011243 

.  132560 

Q 

004616 

.008425 

009587 

.  016022 

-.047575 

.  698254 

M' 

-.033141 

.  025961 

020586 

.  024431 

-.  049677 

.425817 

N' 

-.017248 

, 038471 

-.051417 

.  085977 

225999 

2. 427820 

P' 

-.  005173 

. 003631 

003042 

.003877 

008676 

.  115150 

Q' 

-.001135 

.  004415 

006753 

. 012502 

-.036708 

.  606979 
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TABLE  VII— STATION  NUMBERS 
n  =  7 


1 

2 

3 

4 

5 

6 

7 

M 

1 

1. 243487 

0. 376396 

-0. 026266 

0. 009112 

-0. 005896 

0. 006025 

-0. 006513 

N 

.  667840 

. 480602 

-.063889 

.  031590 

-.  026481 

. 033195 

-.042160 

P 

.218415 

.  022882 

003069 

.001211 

-.000353 

.000925 

-.000363 

0 

.053049 

.  027042 

-.006769 

. 003890 

-.003599 

.004829 

-.  005357 

M' 

.  702228 

.458303 

-.050122 

.  019989 

-.012820 

. 011370 

-.011099 

A' 

. 454474 

.  597757 

-.  122429 

.  069352 

-,  057533 

. 062529 

-.  072080 

P' 

. 164382 

.  032358 

-.  006038 

.  002679 

-.001831 

.  001689 

-.001614 

Q' 

.042465 

.038456 

-.013441 

.003611 

007724 

.  003815 

-.  010037 

M 

2 

-1.321299 

0.  905437 

0. 446479 

-0. 055674 

0. 029812 

-0. 027896 

0.  028701 

N 

-.  570270 

1.  435730 

1.028397 

-. 192270 

. 133730 

-.  153603 

.  185730 

P 

-.154452 

. 182773 

. 031489 

-.007052 

.004233 

-.004239 

.003780 

Q 

-.031847 

.  228718 

.  068932 

-.022639 

.017853 

-.022136 

.023463 

M' 

-.357636 

. 764868 

.  485193 

-.071706 

.  038132 

031029 

. 028988 

N' 

191650 

1. 234950 

1. 123273 

247832 

.  170920 

-.170555 

. 188216 

P' 

-. 058807 

.166641 

.036328 

-.009168 

.005346 

-.004565 

. 004199 

O' 

-.013146 

.209297 

.  079621 

-.029438 

.022542 

-.023823 

.026110 

M 

3 

1.672922 

-0.  511106 

0.  737737 

0.  516672 

-0. 104856 

0. 081487 

-0. 077078 

N 

.  701415 

-.  751768 

1.931045 

1.  718603 

-. 468955 

.448232 

-.  498585 

P 

.  185835 

-.  069049 

. 162183 

.040990 

-.014225 

. 012162 

-.  010056 

0 

.  037666 

-.083877 

. 366024 

.  130958 

-.069956 

.063490 

-.062408 

M' 

.  355047 

-.329159 

.  692136 

.532099 

108454 

.073511 

-.  063669 

N' 

.186093 

-.492459 

1.819566 

1.  771836 

484684 

. 403671 

-. 413268 

P' 

. 056122 

048434 

.  156613 

.  042913 

-.014528 

. 010624 

-. 009166 

O' 

.012374 

-.  059075 

. 353730 

.  137128 

-.061218 

.  055439 

056987 

M 

4 

-1.605312 

0-  409927 

-0.  250374 

0. 629063 

0.  592219 

-0. 182666 

0. 144688 

AT 

-.664780 

. 597643 

625274 

2. 291470 

2.  574726 

-1.002357 

.935220 

P 

-.  174510 

. 052757 

-.037056 

. 147488 

.  051989 

-.026098 

. 018551 

0 

-.035122 

.063907 

-.082279 

.  480975 

.  218352 

136173 

.115110 

M' 

-.317567 

. 247432 

-.216700 

.  623577 

.  584350 

-.  157315 

.115706 

AT' 

-. 164912 

.366953 

-.543417 

2.  273022 

2.  538692 

-.  861833 

.  750618 

P' 

-.049381 

. 034761 

-.033166 

. 147041 

.050504 

—.021773 

.  016464 

O' 

-.010326 

.042284 

-.073703 

. 479557 

. 212061 

-. 113560 

.  102355 

M 

5 

1. 129029 

-0. 264373 

0. 134585 

-0. 136596 

0.  551252 

0. 668960 

-0. 220971 

AT 

.464325 

384003 

. 334325 

-.  480675 

2.  570992 

3. 589325 

-1.425676 

P 

. 121269 

-.033334 

.  019010 

-.021700 

.  136202 

. 063551 

-.027154 

0 

.024312 

-.040333 

.  042147 

-.  069981 

.  580855 

. 330701 

168426 

M' 

. 234133 

-. 168109 

.122950 

-. 143020 

. 568583 

.  634247 

-. 197770 

N' 

. 120971 

-.248408 

.  306705 

-,  503685 

2. 649577 

'  3.397336 

-1.281188 

P' 

. 036034 

-.023166 

.  017933 

-.022914 

. 139031 

.  057851 

027322 

O' 

. 007887 

-. 028149 

.039907 

-.  073907 

. 593022 

.  300912 

169318 

M 

6 

-0. 617807 

0. 137401 

-0. 064103 

0.  054068 

-0.  084474 

0.  507772 

0. 632193 

N 

-. 252961 

. 199169 

-.  158887 

. 189539 

-.383331 

2. 883613 

4. 007039 

P 

-. 065852 

. 017129 

-.  003878 

.003211 

-.  014243 

. 130010 

. 051386 

Q 

-.013170 

. 020712 

019672 

. 026453 

-.  060230 

. 684756 

. 318021 

M' 

-.  125078 

.  086131 

-.  058376 

.  057414 

092005 

. 516450 

. 704705 

AT' 

-.064447 

. 127049 

-.  145336 

.  201491 

-.417397 

2.  931074 

4. 491463 

P' 

-.019183 

.  011759 

003393 

.  008826 

-.  015450 

. 131144 

.  069349 

O' 

-. 004186 

. 014281 

-.  018616 

.  028441 

-.065330 

.  690659 

. 430362 

M 

7 

0.  205449 

-0. 044613 

0. 020059 

-0. 015852 

0. 021382 

-0. 053078 

0.498306 

N 

.083943 

-.064603 

. 049676 

-.  055505 

.  096798 

-. 295079 

3. 334065 

P 

. 021819 

-.005533 

.  002756 

-.  002373 

.  003481 

-.  009551 

. 130932 

0 

. 004358 

-.006689 

. 006106 

-. 007644 

.014713 

-. 049981 

. 820708 

M' 

. 033023 

-.022307 

.  014642 

-.013579 

.  018939 

-.044219 

.  420133 

N' 

.016993 

-.032878 

.  036425 

-.047602 

.  085721 

-.  245636 

2. 816717 

P' 

.005053 

-.  003033 

.  002091 

-.002060 

.  003074 

007854 

.114318 

O' 

.  001102 

003682 

. 004637 

-.  006638 

. 012992 

-.041094 

. 716958 
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'EABLE  VIII— STATION  NUMBERS 


n  =  8 


k 

3 

1 

2 

3 

4 

5 

6 

7 

8 

M 

1 

1.312192 

0. 364019 

-0. 021829 

0. 006490 

-0. 003545 

0. 003081 

-0. 003931 

0. 005039 

N 

. 695399 

.  462793 

-.  052953 

.022453 

-.015893 

.  016958 

-.025621 

.  037691 

P 

.226483 

.  021401 

-.002485 

.000839 

-.000499 

.000464 

-. 000622 

.000684 

Q 

.054447 

.  025255 

-.  005476 

. 002694 

-.002103 

.002422 

-.003869 

.004929 

M' 

.727233 

.  444367 

-.043004 

.  015242 

-.  008698 

. 007082 

-.  007622 

.008334 

N' 

.  467152 

.  577362 

-.  104816 

.  052799 

-.038988 

.038933 

-. 048943 

.062492 

p' 

.168111 

.  030533 

-.005119 

.002003 

-.001220 

. 001045 

-.001145 

.001228 

Q' 

.  043271 

.  036245 

-.011294 

.006433 

-.  005146 

.005456 

-.007123 

.008863 

M 

2 

-1.600390 

0. 955663 

0. 428499 

-0. 045077 

0. 020351 

-0.  016189 

0. 019610 

-0. 024337 

N 

-.682210 

1. 507998 

. 984086 

-.  155335 

.091123 

-.089039 

.127790 

182002 

P 

-.  183160 

. 188781 

.029123 

005560 

.002808 

002412 

.003082 

-.003289 

Q 

-.037524 

.  235968 

.  063696 

-.017807 

.011838 

-.  012593 

. 019183 

-.023705 

M' 

-.415706 

. 797175 

. 468735 

-.060778 

.  028715 

021414 

.  021627 

-.023284 

N' 

-.221089 

1.282229 

1.082551 

-.209723 

. 128559 

-.117660 

.140686 

174660 

P' 

-.  067464 

.170868 

. 034087 

-.007611 

.003953 

-.003129 

.003272 

-.003421 

Q' 

-.015018 

.214419 

. 074659 

-.024428 

. 016664 

-.016330 

.020358 

-.024694 

M 

3 

2.  337500 

-0.  630523 

0. 780389 

0. 491625 

-0. 082648 

0.  054514 

-0.060500 

0. 071477 

N 

.  967933 

-.923581 

2.036155 

1. 631308 

368956 

.299530 

-.394067 

.  534419 

P 

.  254178 

-.083331 

. 167794 

. 037442 

010886 

.  007975 

-.009418 

,  009602 

Q 

.051181 

-. 101109 

. 378439 

.  119546 

-.045861 

.  041626 

-.058605 

.069206 

M' 

.462586 

-.388839 

. 722423 

.512112 

-.091426 

.056607 

-.  052469 

.054058 

N' 

.240602 

-.  579786 

1. 894495 

1.  702145 

-.408100 

.  310711 

-.341167 

.  406216 

P' 

. 072148 

-.056240 

.  160734 

. 040069 

-.012016 

.008120 

-.  007864 

.007911 

Q' 

.015838 

-.068533 

.362854 

. 127977 

-.050617 

.  042369 

-.048918 

. 067107 

M 

4 

-2. 630070 

0.  593584 

-0. 315718 

0. 667195 

0.  558819 

-0. 143453 

0. 132430 

-0. 143916 

N 

-1.075644 

.861856 

-.786292 

2.  424356 

2.  424356 

-.786292 

.861856 

-1.075644 

P 

-.279850 

.074709 

-. 045646 

. 152881 

.  046981 

-.020073 

.  020253 

-.019152 

Q 

-.  055949 

. 090392 

-.101283 

.  498324 

. 197214 

-.  104712 

.  126011 

-. 138032 

M' 

-.466670 

.329868 

-.258293 

.  650770 

.  561588 

135711 

.  104892 

-.  100193 

-.240468 

. 487559 

-.646305 

2. 367821 

2. 436346 

-.743106 

. 681473 

-.750848 

P' 

-.071591 

. 045534 

-.038818 

. 150902 

. 047159 

-.018614 

. 015449 

-.014572 

<?' 

-.015627 

. 055337 

-.086223 

.  491981 

. 197950 

-.097075 

.096084 

-.105178 

M 

5 

2. 192995 

-0.  454018 

0.  201514 

-0. 175130 

0. 584107 

0. 633389 

-0. 237005 

0. 215982 

N 

. 890699 

-.656778 

. 499217 

-.614926 

2.  718847 

3. 393592 

-1.  539301 

1. 613222 

P 

. 230544 

-.055975 

. 027793 

-. 027134 

. 141097 

.058226 

-.034728 

.028253 

Q 

. 045912 

-.  067648 

. 061577 

-.087460 

.601511 

.302906 

-.215987 

.203600 

M' 

. 389428 

-.253524 

.  165704 

-.  170613 

. 591124 

.  614158 

-.191506 

. 155413 

N' 

. 199636 

-.  373346 

. 412443 

-. 599857 

2.  750891 

3.  287044 

-1.241724 

1. 164103 

P' 

. 059201 

-.034316 

. 023783 

-. 026822 

. 142374 

.054969 

-.027027 

.022340 

Q' 

.  012882 

-.041657 

. 052749 

086477 

. 606915 

.285874 

-.168026 

. 161236 

M 

6 

-1.350286 

0.  265559 

-0. 108064 

0. 078152 

-0. 102878 

0.  520374 

0. 709714 

-0. 274441 

N 

-.546005 

. 383370 

-.267117 

.  273370 

-.  466005 

2. 952258 

4. 523995 

-2. 046630 

P 

-. 140865 

. 032367 

-.014611 

.011568 

-.016907 

.  131625 

.  070024 

-.034431 

Q 

-.027983 

.  039093 

-.032351 

. 037251 

-.071468 

.692634 

.434640 

-.248035 

U' 

-.263362 

. 163432 

-.098083 

. 084195 

-.115917 

.  543778 

.662638 

-.230699 

N' 

-. 134575 

.  240192 

-.  243600 

. 294907 

-.525024 

3. 081687 

4. 216400 

-1.725908 

P' 

-.039810 

.021885 

-.013810 

. 012656 

-.019020 

. 135367 

. 062213 

-.032195 

Q' 

-.008646 

. 026552 

-.030611 

.  040765 

-.080395 

. 712714 

.385904 

-.232314 

M 

7 

0. 654484 

-0. 124395 

0. 048081 

-0. 031856 

0.034801 

-0. 066829 

0. 484484 

0. 662747 

N 

.  263837 

179348 

. 118705 

-. 111252 

. 157170 

-.  370670 

3. 239551 

4. 867795 

P 

.  067914 

-.015052 

. 006430 

-.004626 

. 005472 

-.011570 

.126311 

.055545 

Q 

.  013468 

-.018172 

.  014232 

014889 

.023114 

-.060517 

.791745 

.  399381 

M' 

.  130945 

-.079057 

. 045388 

036054 

.041862 

-.077799 

.  500211 

.  723627 

N' 

.  066783 

-.  116069 

.112606 

-.  126105 

. 189086 

431248 

3.341811 

5. 337321 

P' 

. 019727 

-.010524 

. 006331 

-.005328 

. 006593 

-.013324 

.128649 

.072134 

Q’ 

.004280 

-.012763 

. 014029 

-.017155 

.027846 

-.069680 

. 806281 

. 519815 

M 

8 

-0.  202414 

0. 037821 

-0.014268 

0. 009104 

-0.009307 

0.  015388 

-0.045167 

0. 487926 

N 

-.081470 

.064494 

-.035205 

.031774 

-.041992 

.085181 

296574 

3.  754721 

P 

-.020947 

.004560 

-.001898 

.001312 

-.001443 

. 002573 

-.008294 

. 129520 

Q 

-.004150 

. 005504 

-.004201 

. 004221 

-.006093 

.  013452 

-.051711 

. 941457 

M' 

-.033107 

. 019722 

-.011093 

. 008527 

-.009326 

. 015058 

-.039785 

. 414996 

N' 

-.016869 

.028937 

-.027507 

.029809 

-.042089 

.083310 

-  260993 

3. 198326 

P' 

-.004979 

. 002618 

-.001541 

. 001262 

-.001450 

.  002493 

-.007173 

.113558 

«' 

-.001080 

. 003174 

-.003414 

.004029 

-.006122 

. 013034 

-.044710 

.  825790 
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TABLE  IX— COMPARISON  OF  RESIiLTS 


Number 

of 

stations 

Torsion 

Bending 

Coupled 

V  Co 

\  Co 

’V  Co 

V  ^0 

V  Bo 

V  Bo 

V  Bo 

^2-*  - - 

V  Bo 

/mo? 
^3-1/ - 

V  .Bo 

Station-Function  method 

1 

2 

3 

1.549 

1.571 

1.571 

3. 493 

3. 516 

3. 516 

3. 46 

3. 48 

3.  48 

4. 626 

4. 689 

21.71 

22. 04 

19.7 

20.6 

7.  502 

60. 20 

48.2 

Weighted  influence  coefflcients 

2 

4 

1.675 

1.671 

5.39 

4.73 

3.56 

3. 52 

15. 63 

22.80 

Exact  theoretical  value 

1.571 

4.  712 

7.  854 

3. 516 

22. 04 

61.70 

3.49 

20.6 

49.1 

TABLE  X— STATIONS  REQUIRED  FOR  SATISFACTORY  ACCURACY 


Method 

Torsion 

Bending 

V  Co 

"Vl; 

V  Bo 

V  Bo 

V  -Bo 

Station  Functions . 

Weighted  influence  coefficients _ _ 

1 

2 

3 

4 

4 

1 

3 

2 

6 

3 

M.  S.  GOVERNMENT  PRINTING  OFFICE:  i»St 


Axis 


Moment  about  axis 


Angle 


Velocities 


Force 
(parallel 
to  axis) 
symbol 

_ _j 

Designation 

Sym¬ 

bol 

Designation 

Sym¬ 

bol 

Positive 

direction 

Designa¬ 

tion 

Sym¬ 

bol 

Linear 
(compo¬ 
nent  along 
axis) 

Angular 

LongitiiHinnl 

X 

X 

L 

Y - >Z 

Roll 

e 

u 

P 

T.fttersl 

Y 

Y 

M 

Z - >X 

Pit.cli 

V 

q 

Normal 

z 

z 

N 

X - >Y 

4' 

w 

r 

Absolute  coefficients  of  moment 

p _ p  _ ^ 

^‘~qbS  '"’^-qcS 

(rolling)  (pitching) 


^’'~qbS 

(yawing) 


Angle  of  set  of  control  surface  (relative  to  neutral 
position),  5.  (Indicate  surface  by  proper  subscript.) 


4.  PROPELLER  SYMBOLS 


D 

P 

PID 

V' 

V, 

T 

Q 


Diameter 
Geometric  pitch 
Pitch  ratio 
Inflow  velocity 
Slipstream  velocity 


Thrust,  absolute  coefficient  Ct— 
Torque,  absolute  coefficient  Cq= 


T 

pnW* 

pn^I^ 


P 


C. 


V 

n 


Power,  absolute  coefficient  Cp= 


P 


Speed-power  coefficient = 


Efficiency 

Revolutions  per  second,  rps 


Effective  helix  angle = tan  ‘ 


5.  NUMERICAL  RELATIONS 

1  hp=76.04  kg-m/s=’550  ft-lb/sec  1  lb=0.4536  kg 

1  metric  horsepower=0.9863  hp  1  kg=2.2046  lb 

1  mph=0.4470  mps  1  mi=  1,609.35  m=5,280  ft 

1  mps=2.2369  mph  1  m=3.2808  ft 


